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Recent work on random field Ising model is described briefly emphasising exact solutions of 
the model in simple cases and their relevance in understanding equilibrium and non-equilibrium 
properties of systems with quenched disorder. 



I. INTRODUCTION: 

Extended systems with quenched random disorder of- 
ten possess a large number of nearly degenerate states 
separated by high energy barriers. This kind of energy 
landscape gives rise to very complex relaxation phenom- 
ena at low temperatures as well as hysteresis in the sys- 
tem. Even a weak disorder may modify the system sig- 
nificantly and destroy the long range order in the equi- 
librium state of the system. One needs a simple theo- 
retical model to make sense of the rich, complex, and 
vast amount of experimental data in this field. The ran- 
dom field Ising model is perhaps the simplest minimal 
model that fits the bill. Although very simple to state, 
the model is not exactly solvable except in a few special 
cases. Computer simulations of the model are helpful, 
but are often unable to resolve the key questions because 
they suffer from metastability and slow relaxation rates 
in the system just as much as the laboratory experiments. 
For specificity, we focus on disordered magnets. Some of 
the important questions one would like to be answered 
are the following: What is the lower critical dimensional- 
ity below which any disorder, no matter how small, would 
destroy the long range order in the system at zero tem- 
perature? Could the disorder cause a first order jump 
discontinuity in the response of the system to an applied 
field ? If so, in what circumstances? In recent years, 
we have examined these questions in some exactly solved 
cases of the model. Due to the limitation of space, we 
only indicate the scope and status of our studies, and 
mention the results obtained so far. Interested readers 
should consult the references for more technical details. 



II. THE MODEL 

The model is defined on a lattice. Each site is labeled 
by an integer i, and carries an Ising spin Si ( Si = ±1 ), 
a quenched random magnetic field hi, and an externally 
applied uniform field h. The quenched fields {hi} are in- 
dependent identically distributed random variables with 
a continuous probability distribution (f>(hi). For conve- 
nience, we assume (f>(hi) to be Gaussian with mean value 
zero, and variance a 2 . The nearest neighbor interaction 
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can be ferromagnetic (J>0) or anti-ferromagnetic (J<0). 
The Hamiltonian of the system may be written as, 

i,j i i 

The spin Si experiences a net field fi on it that is given 
by, 

f i = jJ2S j + h i + h (2) 

3 

The Glauber dynamics of the system at temperature 
T is specified by the rate Ri at which a spin Si flips to 
-Si 

Ri = - [1 - S t tanh{/ 4 /(fc B T)}] , (3) 

T 

Here r sets the basic time scale for the relaxation of 
individual spins. The energy of the spin Si is equal to 
—fiSi. If fi and Si have the same sign, we say that the 
spin is aligned along the net field at its site. The energy 
of a spin is the lowest if it is aligned along the net field 
at its sight. We are interested in the properties of the 
model at zero temperature ( T = ), and on time scales 
much larger than t. In this limit, the dynamics simplifies 
to the following rule: choose a spin at random, and flip 
it only if it is not aligned along the net field at its site. 
Repeat the process till all spins are aligned along the net 
fields at their respective sites. The dynamics described 
above always brings the system to a fixed point state that 
is stable against single spin flips. The fixed point corre- 
sponds to a local minimum of the energy of the system. 
The system possesses a thermodynamically large number 
of fixed point states each marking a domain of attraction 
in the phase space of the system. The particular local 
minima that the system reaches depends on the domain 
that contained the initial state of the system. 

Although the fixed points are stable, they can be 
thought to represent the metastable states of the sys- 
tem in the sense that they are local minima of energy. 
In this context, it is useful to consider the time scales of 
interest in physical systems, and the appropriateness of 
the zero-temperature dynamics as a model. In a typi- 
cal experiment on non-equilibrium behavior of a physical 
system, there are at least four time scales: (i) time r that 
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an individual spin takes to relax, (ii) time n that the sys- 
tem takes to relax to a metastable state, (iii) time r 2 over 
which the applied field changes, and (iv) life time T3 of 
the metastable state. In complex systems, t\, t%, and T3 
may each contain an entire spectrum of time scales. In 
physical systems relevant to our model, the shortest time 
is r which may be taken to be unity to set the scale. The 
next larger time is T\ — v x r, where v is the number 
of iterations of the dynamics to reach a fixed point. The 
applied field is assumed to vary very slowly (driving fre- 
quency goes to zero!) so that r 2 >> r%. Specifically, it 
means that the applied field is held constant during the 
relaxation of the system. The time T3 is infinite. Thus, 
the present model is applicable if r 3 >> t 2 >> n >> r. 
These conditions fit a wide class of complex magnetic 
materials that have a large number of metastable states 
separated from each other by barriers much larger than 
the available thermal energy. 



III. RELEVANCE TO EXPERIMENTS 

At an intuitive level, quenched random fields may be 
thought to arise from impurities and imperfections in dis- 
ordered magnetic materials. A more precise connection 
between random field Ising model and experiments was 
made by a theoretical argument due to Fishman and 
Aharony Q suggesting that the critical behavior of a 
weakly diluted anti-ferromagnet in a uniform external 
field should be in the same universality class as that 
of a ferromagnet in a random external field. The es- 
sential idea is that an anti-ferromagnet without dilution 
has two sub-lattices on which the spins are oriented op- 
posite to each other. One of the sub-lattice is aligned 
with the applied field. In the presence of dilution, lo- 
cally the sub-lattice with most spins tends to align with 
the applied field in competition with the global anti- 
ferromagnetic order in the absence of dilution. As a re- 
sult the applied field acts as an effective random field 
coupling to the anti-ferromagnetic order parameter i.e. 
the staggered magnetization. The effective random field 
produced by the applied field is proportional to the ap- 
plied field. The strength of the random field is there- 
fore easily controlled. One can do scaling studies with 
varying strengths of disorder in a system by simply tun- 
ing the applied field rather than making fresh samples 
with different degrees of dilution. This feature has a 
far reaching significance in investigating the theoretical 
model experimentally, and the interaction between the- 
ory and experiment has helped the field develop con- 
siderably, and clarified many fine points of the model 
0, RbiCoF^ is a good two dimensional Ising anti- 
ferromagnet. It consists of layers of magnetic ions with 
a single dominant intralayer exchange interaction, and 
an interlayer interaction which is smaller by several or- 
ders of magnitude. It is very anisotropic so that the 
spins can be well represented as Ising spins. The mate- 
rial can be magnetically diluted by introducing a small 



fraction of manganese ions in place of cobalt ions. Thus 
crystals of RbiCo x Mn\_ x F± are good examples of a two 
dimensional diluted anti-ferromagnet suitable for an ex- 
perimental realization of the two dimensional random 
field Ising model. In three dimensions, the most stud- 
ied dilute anti-ferromagnet is Fe x Zri\- x Fi crystal. In 
the pure ferrous fluoride FeF2 crystal, the ferrous ions ( 
Fe ++ ) are situated approximately on a body centered 
tetragonal lattice. Each ferrous ion is surrounded by a 
distorted octahedron of fiurinc ( F~ ) ions. The predom- 
inant interactions are a large single-ion anisotropy, and 
an anti-ferromagnetic exchange between nearest neigh- 
bor ferrous ions. The magnetic moments of ferrous ions 
on the corners of the tetragonal cell are anti-parallel to 
the magnetic moments of Fe ++ ions on the body cen- 
ters. The large crystal field anisotropy persists as the 
magnetic spins are diluted with Zn. The diluted crys- 
tal remains an excellent Ising anti-ferromagnet for all 
ranges of magnetic concentration x. Furthermore crys- 
tals with excellent structural quality can be grown for 
all concentrations x with extremely small concentration 
variation 5x < 10~ 3 . These attributes combine to make 
Fe x Zn\- x Fi the popular choice for experiments on di- 
luted anti-ferromagnets, although experiments have been 
done on several other materials as well H. 



IV. DOMAINS IN RANDOM FIELDS 

In one of the earliest studies of the random field Ising 
model, Imry and Ma |l( argued that quenched random 
fields in a system may cause a uniform ferromagnetic 
state to break into domains. The argument is essentially 
as follows. Consider a uniform ferromagnetic state at 
zero temperature. In the presence of random fields, a 
strategically placed domain of linear size L may turn over 
and gain an energy of the order of uL d / 2 .. However, this 
would create a domain wall that would cost an energy 
of the order of 3L d - x . If d/2 > (d - 1), i.e. if d < 
2, then for any a, there will be a characteristic length 
over which the bulk energy gain will overcome the cost 
of the surface energy. In other words, domains will occur 
spontaneously if d < 2. The above argument is intuitively 
appealing, but nonrigorous. It is also inconclusive at the 
lower critical dimensionality because the gain and the 
cost of energy both scale linearly with L at d = 2 . The 
situation in two dimensions can be clarified by taking into 
account the roughness of the domain wall. The work of 
Binder p| supplemented with numerical simulation of 
a toy model predicts that the gain in energy scales as 
uLlogL in d = 2. Thus the domain argument predicts 
the absence of long range order in two dimensions. 

It remained unclear for several years if the results ob- 
tained from the domain argument were correct. The con- 
troversy was generated by the existence of another argu- 
ment based on a field theoretic method (dimensional re- 
duction) that predicted the lower critical dimensionality 
of random field Ising model to be three. Initially, com- 
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puter simulations of the model as well as experimental 
observations did not help to resolve the controversy. It 
took several years to realize the error in the dimensional 
reduction argument. The difficulty with both theory and 
experiment was in the interpretation of results. Concen- 
tration gradients in the diluted anti-ferromagnetic sam- 
ples tend to round off a transition and affect the mea- 
surements of critical behavior drastically. Further, the 
majority of experiments were performed on samples pre- 
pared in two separate ways; (i) cooling the sample in zero 
magnetic field, and (ii) cooling it in a magnetic field, and 
turning the field off at the end. Experiments on samples 
prepared in the two ways yield different results. Three 
dimensional field cooled samples show no long range or- 
der and were first thought to show that the lower critical 
dimensionality of the model is three. But three dimen- 
sional samples cooled in zero field showed long range or- 
der. It was only after several years of controversy that 
the experimental situation resolved itself in favor of the 
domain argument, i.e. no long range order in two di- 
mensions, but long range order in three dimensions. The 
main point that was clarified by theoretical work is that 
the field cooled state is not an equilibrium state. It re- 
laxes logarithmically slowly, and one should not expect to 
see an equilibrium ordered state in a field cooled sample 
over any reasonable experimental time scale. 



V. HYSTERESIS 

Permanent magnets are typically a two phase solid ma- 
terial with fine magnetic particles of one phase embedded 
in the other phase. The precipitation is carried out in a 
magnetic field, and needle like magnetic particles are ori- 
ented with their long axis parallel to the field direction. 
If the field is reversed after the material is set in a solid 
matrix, the needles can not physically turn over as in 
a liquid matrix. Instead, the magnetic domains inside 
the particles have to reverse themselves. This requires a 
threshold field that varies from particle to particle. Ex- 
perimentally observed magnetization of the material in a 
smoothly increasing applied field may look smooth on a 
macroscopic scale but on a microscopic scale, it is made 
of steps of irregular widths and heights. This is known as 
Barkhausen noise. Hysteresis in such materials is some- 
what different in character from the usual hysteresis that 
arises from a delayed response of the system to the driv- 
ing field. In the present case the hysteresis is dominated 
by the energy landscape of the system. There are a very 
large number of local minima with barriers of varying 
heights between them. Suppose the applied field is in- 
creased sufficiently so that the system jmps over a barrier 
from one minimum A to another minimum B. If the field 
is then reversed, it may be more favorable for the sys- 
tem to jump over a smaller barrier to another minimum 
C rather than return to A. This is the kind of hystere- 
sis that is modeled very well by the zero temperature 
dynamics of the random field Ising model Q ■ We have 



solved the ferromagnetic model exactly on a Bethe lattice 
of an arbitrary coordination number z, and obtained ma- 
jor as well as minor hysteresis loops, and power laws for 
the Barkhausen noise analytically [g, 0, Hoi ITTl IT^L H"^ . 
In the following, we describe our main results without 
going into the details of calculation. 

The main limitation of the method used to obtain the- 
oretical hysteresis loops has been its restriction to Bethe 
lattices, and to initial states having saturated magneti- 
zation i.e. all spins pointing in the same direction. More 
recently the method has been extented to the case of a 
random initial state [l6| but we shall not go into it here. 
Suppose all the spins are pointing down. This is a stable 
state in an applied field h = — oo. We need to calculate 
the magnetization as the applied field in increased from 
h = — oo to ft, = oo. The other half of the hysteresis loop 
can be obtained by a symmetry of the Ising model, and 
that of the Gaussian distribution centered at the origin. 
In the first instance, we calculated the magnetization on 
a one dimensional lattice @- The calculation utilized 
the fact that the evolution on a linear segment of the lat- 
tice contained between two sites that flip up before their 
nearest neighbors is independent of the rest of the lat- 
tice. It was not immediately obvious how the calculation 
could be generalized to a Bethe lattice. What makes the 
generalization possible is the following: Suppose the spin 
Si in the deep interior of the Bethe lattice does not flip 
up as the applied field increases from ft, — oo to ft,, then the 
evolution on each of the z sublattices that meet at site 
i is independent of each other upto the applied field ft. 
The second point is that the stable state at applied field 
ft, does not depend on the order in which the spins are re- 
laxed if the interactions are ferromagnetic. This property 
of the ferromagnetic model is called the Abelian property. 
It allows us to relax the nearest neighbors of site i inde- 
pendently of each other and before site i is relaxed. The 
method involves calculating the conditional probability 
that a nearest neighbor of site Si is up before Si flips 
up. This conditional probability is obtained as a fixed 
point of a recursion relation describing the relaxation of 
sites starting from the boundary of the Bethe lattice and 
working in. 

Hysteresis in the ferromagnetic model on a Bethe lat- 
tice has symmetrically placed jump discontinuities on the 
two halves of the hysteresis loop if a is smaller than a 
critical value er c , and the coordination number z of the 
lattice is greater than or equal to four. In other words, 
hysteresis loops on a Bethe lattice of coordination z > 
3 are qualitatively different from those on lattices with 
z=2, and z=3. For z < 3 there is no jump discontinuity 
in the hysteresis loops for any amount of disorder. For 
z > 4, there is a critical value of a that characterizes the 
Gaussian random field distribution. If a is less than the 
critical value a c , the magnetization in increasing field has 
a macroscopic first-order jump at an applied field h c > 
J. As a increases to er c , h c decreases to J, and the first- 
order jump in magnetization reduces to zero. The system 
shows non-equilibrium critical behavior at ft = h c , and 
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a = a c . For z=4, cr c = 1.78 approximately. The value of 
critical disorder a c increases with the coordination num- 
ber z of the lattice. At the critical point the Barkhausen 
noise shows a power law distribution. The probability of 
avalanches of size s scales as s~ 3 / 2 . Nonequilibrium crit- 
ical point phenomena and its relationship with the coor- 
dination number of the lattice appears to hold on other 
lattices as well. Numerical simulations and theoretical 
arguments on several periodic lattices embedded in two 
and three-dimensional space show that hysteresis on pe- 
riodic lattices with z > 4 is qualitatively different from 
that on lattices with z < 4. Although there are some 
similarities between Bethe lattices and periodic lattices 
of the same coordination number, there are differences as 
well. The differences are related to the bootstrap perco- 
lation instability on some periodic lattices |15| . 

Minor hysteresis loops have also been obtained, and 
their calculation reveals an interesting feature of the 
model that has an experimental significance. Consider 
two halves of the major hysteresis loop connecting states 
of saturated magnetization in opposite directions. If the 
applied field is reversed while the system is on one of 
these halves, the magnetization trajectory branches off 
and heads towards the other half of the major loop. We 
find that the trajectory in reversed field meets the other 
half of the major loop exactly when the field has been 
reversed by an amount 2J irrespective of the point of 
reversal. This result provides an interesting possibility 



for measuring the exchange interaction J in a hysteresis 
experiment. 

We have also studied hysteresis in the anti- 
ferromagnetic random field Ising model at zero temper- 
ature |14| . Hysteresis in the anti- ferromagnetic model 
is qualitatively different from that in the ferromagnetic 
model because it does not show Barkhausen noise. On ac- 
count of the anti-ferromagnetic interactions, a spin turn- 
ing up in an increasing field blocks its nearest neigh- 
bors from turning up. Thus there is no microscopic 
avalanche of up-turned spins as in the case of ferromag- 
netic interactions. A spin that turns up in increasing field 
in the anti-ferromagnetic model, occasionally causes its 
nearest neighbor (that had turned up earlier) to turn 
down. As the applied field increases from — oo to +oo, 
a small fractions of sites flip three times, first up, then 
down, and finally up again. Also, the dynamics of the 
anti-ferromagnetic model is non-Abelian. These features 
make the anti-ferromagnetic dynamics rather complex, 
and an exact solution of the model becomes difficult. So 
far the anti-ferromagnetic model had been solved exactly 
in one dimension only, and that too for a rectangular dis- 
tribution of the random field of width 2 A, where A < \J\. 
Recently, we have been able to extend the calculation of 
hysteresis in a one dimensional random field Ising model 
to an arbitrary continuous distribution of the random 
field 0. 
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